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Abstract 


The present study was conducted to investigate the trend of land use changes and its impact on 
erosion and sediment, using Landsat satellite images of 1987, 2010 and 2019 and EPM model in 
the Kal-Aji watershed of Golestan province. Land use map was prepared using the maximum 
likelihood method. Training data were used to evaluate the accuracy of the results. Then, the most 
important methods of accuracy assessment, including total accuracy and classification kappa 
coefficient, were extracted. Those for 1987, 2010 and 2019, the total accuracy were 95, 98 and 
95 and kappa coefficient were 0.93, 0.96 and 0.93, respectively. Erosion and sedimentation zoning 
maps were performed using the resulted land use maps, and such factors as slope, lithology, soil 
and precipitation were performed using the EPM method. The results showed that one of the most 
important factors in the Kal-Aji watershed is the loss of rangeland vegetation and its conversion 
into agricultural land. According to the hypothesis, the role of land use changes in the amount and 
increase of erosion and sediment was significantly determined. That is, with the change of land 
use, erosion and sediment have also increased. The reasons for the increase in the amount of 
sediment in the Kal-Aji area are the development of agricultural land uses and the decrease in 
rangeland between 1987 and 2019 (16.83%). These changes have reduced the important role of 
rangelands and vegetation of Kal-Aji watershed and have increased the amount of sediment in the 
outflow of sub-watersheds. Based on the obtained results, the total amount of special sediment 
was related to the period of 1987 (4353 cubic meters), which increased to 5164 cubic meters in 
2019, which indicates the changes in land use that occurred in the watershed. 
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مقاله پژوهشی 


ارزیابی تأثیر تغییرات کاربری اراضی روی فرسایش و رسوب با استفاده از مدل EPM‏ 
(مطالعه موردی: حوضه کال آحی استان گلستان) 
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حکده 


نوع و شدت فرسایش خاک در یک منطقه. تابع شرایط اقلیمی پستی‌وبلندی زمین» خاک و کاربری اراضی 
است که در این cle‏ اهمیت کاربری اراضی به دلیل نقش مؤثر انسان بر OF‏ نسبت به دیگر عوامل زیادتر 
است. هدف پژوهش حاضر بررسی روند تغییرات کاربری اراضی/ پوشش گیاهی و SU‏ آن بر فرسایش و 
رسوب. با استفاده از تصاویر ماهواره‌ای سال‌های AY VI‏ ۱۳۸۹ و ۱۳۹۸ لندست و مدل EPM‏ در حوضه 
کال آجی استان گلستان است. نقشه کاربری اراضی با استفاده از روش حداکثر احتمال تهیه گردید. برای 
ارزیابی صحت نتایج» از داده‌های آموزشی استفاده به عمل آمد. سپس مهم‌ترین روش‌های ارزیابی صحت؛ 
شامل دقت کل و ضریب GUIS‏ طبقه‌بندی استخراج شد که برای سال‌های ۸۱۳۲۳ ۱۳۸۹ و ۱۳۹۸ به ترتیب؛ 
برای دقت کل ٩۵‏ ۹۸ و ۹۵ و ضریب LE‏ ۰۰/۹۳ ۰/۹۱ و ۰/۹۳ درصد بود. نقشه‌های پهنه‌بندی شدت 
فرسایش و رسوب با استفاده از نقشه‌های کاربری اراضی حاصل و عواملی چون شیب. لیتولوژی. خاک و 


بارش با استفاده از روش EPM‏ انجام شد. نتایج نشان می‌دهد از مهم‌ترین عوامل موثر در افزايش فرسایش 
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و رسوب در حوضه آبخیز کال‌آجی از بین رفتن پوشش گیاهی مرتعی و تبدیل OT‏ به اراضی کشاورزی 
است؛ یعنی با تغییر کاربری بالطبع فرسایش و رسوب نیز افزايش يافته است. توسعه سطوح کاربری اراضی 
کشاورزی و کاهش سطوح مرتعی در بین سال‌های ۱۳۹۱ تا ۱۳۹۸ به میزان ۱2/۸۳ درصد نسبت به سال 
۲ است. این تغیبرات نقش مهم تخریب مراتع و پوشش گیاهی حوضه آبخیز کال‌آجی در افزایش میزان 
رسوب خروجی OLE‏ می‌دهد. بر اساس نتایج به دست‌آمده مجموع مقادیر رسوب ویژه مربوط به بازه 
زمانی ۱۳۹۲ EYOF)‏ مترمکعب) بوده که این میزان در سال ۱۳۹۸ به OVIE‏ مترمکعب افزایش داشته است 


که نشان از تغییرات کاربری اراضی رخ داده در حوضه انا 


کلیدواژه‌ها: فرسایش و رسوب» پوشش گیاهی/ کاربری اراضی. GIS EPM‏ کال‌آجی. 


۱- مقدمه 

معضل فرس‌ایش با کاهش حاصلخیزی خاک باعث تخریب اکوسیستم‌های طبیعی نظیر جنگل‌ها و مراتع و 
اکوسیستم‌های کشاورزی می گردد. مطالعات متعددی نشان داده که در میان عوامل مسبب فرسایش خاک کاربری 
اراضی از اهمیت بیشتری نسبت به ساير عوامل برخوردار است. پوشش گیاهی به‌عنوان یک عامل کلیدی در تعیین 
کمیت و کیفیت فرسایش خاک و نیز پاسخ‌های هیدرولوژیکی سطحی زمین است (آدیلیا و همکاران4 ۲۰۱۱ نوع و 
شدت فرسایش عموماً تابع اقلیم. توپوگرافی» خاک و کاربری اراضی است که در این میان SE‏ کاربری اراضی به 
دلیل نقش انسان بر چرخه طبیعت مهم‌تر از سایر عوامل است. نوع کاربری اراضی به همراه شرایط مدیریت 
بهره‌برداری از آن به‌عنوان یکی از عوامل مهم و تأثیرگذار بشری بر روی منابع طبیعی و بخصوص فرسایش خاک به 
شمار می‌آید؛ بنابراین دانستن اينکه هرکدام از کاربری‌های اراضی چه نقشی در تولید فرسایش و رسوب دارند و با 
توجه به روند قبلی و شرایط زمین محیطی منطفه. عوامل مؤثر و محرک در تغییر کاربری اراضی کدام هستند از 
ابتدایی ترین اطلاعات برای مدیران و بهره‌برداران اراضی به شمار می‌آید. با توجه به نقش بهینه کاربری اراضی در 
توسعه پایدار و کاهش هدر رفت منابع طبیعی» می‌توان گفت که تغیبرات کاربری اراضی باعث به هم خوردن چرخه 
bow‏ طبیعی می‌گردد. همچنین به‌تبع آن فرسایش خاک در اثر عدم مدیریت صحیح اراضی و کاربری نادرست آن 
باعث بروز اثرات نامطلوبی بر زندگی اقتصادی- اجتماعی ساکنین حوضه نیز خواهد شد. 

در کشور Ol pl‏ به دلیل عدم توجه به مسئله و قابلیت و تناسب کاربری زمین» بیشتر اراضی به‌صورت نامناسب و 
نامعقول استفاده می‌شوند که این استفاده نادرست به‌شدت فرسایش و رسوب حوضه‌های آبخیز را افزایش می‌دهد. 
مطالعات زیادی در زمینه بررسی تغییر کاربری و نقش آن بر روی فرسایش و رسوب صورت گرفته است که بعضی 


از آن‌ها به شرح زیر است: 


1 Adelia et al 


سال دوازدهم ارزیابی تأثیر تغییرات کاربری اراضی روی فرسایش و .... Wo‏ 


شیلبرگ و همکاران" (THN)‏ به بررسی فرسایش توسط تغییرات استفاده از زمین در کوئینزلند استرالیا پرداختند. 
در این پژوهش میزان فرسایش آبرفتی و ویژگی سیلاب‌های پیش از آن برای ارزیابی تأثیر تخیر کاربری اراضی زمین 
در استرالیا را موردمطالعه قرار گرفته است. در این پژوهش با GPS‏ و عکس‌های هوایی تاریخی در ۱۸ محل در 
سراسر مگافان نرخ فرسایش را به دست آورده و به این نتیجه رسیدند که چرای بیش‌ازحد دام‌ها از پوشش ALS‏ و 
به وجود آمدن رواناب‌ها؛ اختلالاتی را در منطقه به وجود آورده‌اند که موجب فرسایش خاک‌های منطقه شده است. 
پاروسیون و همکاران" (۲۰۱۵) در مطالعه خود با موضوع روشی برای مدل‌سازی اثرات تغییرات آب‌وهوا و کاربری 
اراضی در فرسایش و پایداری خاک در حوضه مدیترانه (لانگدوک» فرانسه) به این نتیجه رسیدند که حساسیت شدت 
فرسایش خاک به تغییرات کاربری اراضی بیشتر از تغییرات اقلیمی است. یوان و همکاران" (۲۰۱۵) در مطالعه خود 
رواناب سطحی و بار رسوب تحت کاربری اراضی مختلف در حوضه کوهستانی تای‌هانگ در شمال چین را 
شبیه‌سازی نمودند. نتایج تحقیق نشان داد که بار رسوب و رواناب در کاربری جنگل و بوته‌زار از همه کاربری‌ها PS‏ 
است. میلوسکی* (۲۰۰۸) خطر فرسایش خاک در پایین‌دست حوضه آبخیز برجالنیکاه در جمهوری مقدونیه را با 
استفاده از تصاویر ماهواره‌ای» روش EPM‏ و ابزار GIS‏ بررسی نمودند و به این نتیجه رسیدند که GIS‏ ابزاری 
قدرتمندی برای پیش‌بینی و تخمین خطر احتمالی فرسایش خاک است. قهرمان‌نژاد و همکارانش (۲۰۱۷) در تحقیقی 
به بررسی تأثیر سناریوهای مختلف مدیریت کاربری اراضی بر میزان فرسایش خاک با مدل USLE‏ در حوضه 
کلیبرچای پرداختند. نتایج کار آن‌ها نشان داد با افزایش اراضی دیم‌زان میزان فرسایش افزایش و با افزایش اراضی 
پوشش درختی» میزان فرسایش کاهش می‌یابد. عابدینی و همکاران (۲۰۱۳) میزان فرسایش و رسوب حوضه آبخیز 
مشکین‌چای (مشکین‌شهر- اردبیل) را با روش EPM‏ برآورد نمودند. نتایج حاصل از پژوهش نشان داد که سه عامل 
مهم توپ وگرافی» لیتولوژی و تغییرات کاربری اراضی در کنترل فرسایش و رسوب نقش مهمی Le!‏ نموده‌اند. مقدار 
ضریب شدت فرسایش dZ)‏ عدد oF‏ میانگین رسوب کل ۱۰۱۲۳۷۹۶ مترمکعب در سال» دبی رسوب ویژه ۱۹۳۹۷/۹ 
و فرسایش ویژه ۲۱۵۵۳/۳ به دست آمد که نشان‌دهنده فرسایش خیلی شدید در منطقه موردمطالعه است. بيات و 
همکاران (۲۰۱۲)» در پژوهشی با استفاده از مدل EPM‏ و به کمک فناوری‌های سنجش زدور و سامانه‌های اطلاعات 
جغرافیایی و استفاده از تصاویر ماهواره *171۷1 به‌منظور تعیین اثر مدیریت کاربری اراضی بر میزان حوضه آبخیز 
دریاچه ارومیه به ارزیابی پرداختند. تاجیک (۲۰۰۷) اثر اقدامات آبخیزداری را بر روی سیل‌خیزی و رسوب‌دهی 
حوضه آبخیز رامیان را بررسی نمود. در بخشی از تحقیق خود با مقایسه نتایج مدل EPM‏ با روش هیدرولوژیکی. 
برآورد بار Glee‏ حد وسط دسته‌های دبی روزانه دریافت» دقت مدل EPM‏ در این حوضه بالغ‌بر ۹۹/۲ درصد است. 

1 Shellberg 
2 Paroissien 


3 Yuan 
4 Milevski 


۳۷۴ جغرافیا و مخاطرات محبطی شمارة سوم 


نبی‌پی لشکریان (۲۰۰۰) به بررسی اثر تغییر کاربری اراضی در میزان فرسایش و رسوب نمود. وی به روش‌های 
و EPM‏ اقدام به مقایسه رسوب‌دهی حوضه در شرایط کاربری اراضی فعلی و آتی (متناسب با استعداد) 
اراضی نمود. نتایج حاصل از این دو روش نشان داد که در صورت اصلاح و احیای پوشش گیاهی و استفاده از اراضی 
بر اساس استعدادشان, فرسایش ویژه در روش MPSIAC‏ به میزان ۹/۲۳ درصد و در روش EPM‏ به میزان ۲۲/۹۸ 
درصد کاهش خواهد یافت. 

اسفندیاری و همکاران (۲۰۱۶) در مطالعه‌ای در حوضه رود ورس قزوین» به بررسی تأثیر تغییر کاربری اراضی بر 
فرسایش پرداختند و به این نتیجه رسیدند که کاربری زراعی از بیشترین و مراتع از کمترین پتانسیل برای فرسایش 
برخوردار بودند. Gale‏ و همکاران (۲۰۱) نیز به ارزیابی تغییر کاربری اراضی مرتعی به دیم‌زار در دو شهرستان رابر 
و ارزوئیه پرداختند. نتایج نشان داد گسترش فعالیت‌های زراعی سبب می‌شود مراتع به زمین‌های کم بازده تبدیل 
شود. فیضی‌زاده (THY)‏ تغیبرات کاربری اراضی و نقش آن در فرسایش حوضه Ob gle‏ را بررسی کردند. آن‌ها 
دریافتند که تبدیل اراضی مرتعی پر تراکم به مراتع کم‌تراکم و زراعت دیم مهم‌ترین عامل در فرسایش منطقه است. 
اصغری‌سراسکانرود و همکاران (۲۰۱۷) نیز به بررسی روند تغییرات کاربری اراضی شهرستان نیر و تأثیر آن بر 
فرسایش منطقه پرداختند. آن‌ها به این نتیجه رسیدند که به ترتیب ۱۲/۵۶ و ۲۷۵۱ درصد از مساحت شهرستان, در 
دو طبقه بسیار پرخطر و پرخطر قرار دارد. به‌طور عمده مناطق با طبقه بسیار پرخطر و پرخطر در هر دو دوره «giles‏ 
در کاربری‌های زراعت دیم و زراعت I‏ - باغات قرار دارند. 

دادرسی‌سبزوار (۱۳۸۶) اثر تغییرات عمق و شیب بر میزان تولید رسوب را موردبررسی قرار داد و نتیجه‌گیری 
کرد که رابطه مستقیمی بین میزان تولید سالانه با زاویه و طول شیب و رابطه عکس با عمق خاک وجود دارد. 
نتایج مطالعات باتنی و گریسمر" (۲۰۰۰) نشان داد افزايش درصد شیب بین ٤‏ الی ۱۱ درصد تأثیر مثبت و 
معنی‌داری بر روی میزان هدررفت خاک دارد. 

حوضه آبخیز کال‌آجی به‌عنوان یکی از زیرحوضه‌های گرگانرود در طی سالیان گذشته فرسایش شدیدی را 
متحمل شده است. این فرسایش به دلایل مختلفی رخ داده که یکی از آن‌ها تغییرات در کاربری زمین و ضعف 
مدیریت منابع آب‌وخاک است. هدف از این مطالعه برآورد تغییرات حاصل شده در فرسایش خاک و رسوب حوضه 
آبخیز کال‌آجی در یک بازه زمانی ۳۲ ساله و بررسی نقش کاربری اراضی و پوشش زمین در میزان فرسایش و رسوب 
با اسفاده از مدل EPM‏ است, برای سنجش تغیرات کاربری زمین از مقایسه تغییرات طبقه‌بندی تصاویر لندست برای 
سال‌های ۰۱۳۱۲ ۱۳۸۹ و ۱۳۹۸ استفاده گردید. 


1 Battanyand Grismer 


سال دوازدهم ارزیابی تأثیر تغییرات کاربری اراضی روی فرسایش و .... Yw‏ 


۲- مواد و روش‌ها 

۱-۲- موقعیت منطقه موردمطالعه 

حوضه آبخیز کال‌اجی با مساحت ۱۷۰۹۹ هکتان ۰/۸۳ درصد مساحت استان و ۱/۶۷ درصد مساحت حوضه 
آبخیز گر گانرود گلستان را به خود اختصاص داده است. این حوضه در محدوده طول شرقی ۵۵ درجه و ۱۸ دقیقه و 
asl ۶‏ تا 00 درجه و ۲۹ دقيقه و ۳۲ انیه و عرض شمالی ۳۷ درجه و ۳۵ adds‏ و acl ٤٤‏ تا ۳۷ درجه و 1۵ دقيقه 
و ۳۷ انیه تقریاً در شمال شهر AIS‏ واقع در شرق استان گلستان است. این حوضه در سرشاخه گرگانرود IB‏ 
داشته و از مناطق سیل‌خیز منطقه است. بلندترین نقطه حوضه VUE‏ متر است که در شمال شرقی حوضه واقع شده 
است و کمترین ارتفاع آن مربوط به خروجی حوضه است که حدود ۲۲۰ متر ارتفاع دارد و بارش متوسط آن 1۲۵ 
میلی‌متر و دمای متوسط سالانه ۱۷۶6 درجه سانتی گراد دارای اقلیم نیمه‌حشک است. بیشترین بارش متوسط کل 
حوضه مربوط به ماه اسفند با 0۰/۱۱ میلی‌متر و کمترین آن به ماه تیر با ۱۲/4۶ میلی‌متر است. شکل شماره ۱ نشان 


دهنده موقعیت محدوده موردمطالعه در تقسیمات سیاسی استان گلستان aal‏ 
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شکل ۱- موقعیت حوضه آبخیز کال‌آجی در تقسیمات سیاسی کشور 


SAPS روش‎ -۲-۲ 


YVA‏ جغرافیا و مخاطرات محبطی شمارة سوم 


تجزیه و پردازش و اعمال تصحیحات لازم بر روی تصاویر ماهواره‌های: ۱- لندست 0 سنجنده TM‏ به تاریخ 
asdi -Y ۵‏ ۷ سنحنده ETM*‏ به تاریخ ۰ -Y‏ لندست ۸ سنجنده OLI‏ به تاریخ ۳۹/۵ 
این روش چهار مشخصه شامل ضریب فرسایش حوضه آبخیز CY)‏ ضریب کاربری اراضی (Ka)‏ ضریب حساسیت 
مقیاس ۰ از سازمان زمین‌شناسی و اکتشافات معدنی کشور تهبه شده است. همچنین» آمار مربوط به 
ایستگاه‌های هواشناسی منطقه در دوره آماری ۳۲ سال آبی ۱۳۹۱ تا ۱۳۹۸ از بایگانی و آمارهای شرکت منابع آب 
ایران و گزارش تلفیق مطالعات منابع آب گرگانرود استخراج گردید. 

—\-Y-Y‏ پردازش تصاویر ماهواره‌ای 

امروزه رایج‌ترین روش برای تهیه نقشه کاربری اراضی و پوشش گیاهی» استفاده از فن‌آوری سنجش‌ازدور است 
(قوشه و همکاران  (Te‏ تصاویر مورداستفاده مربوط به ماهواره لندست و سنجنده‌های OLI ۷۲۲۳ TM‏ 
طبقه‌بندی و استخراج لایه‌های کاربری از تصاویر مربوط به سه دوره انجام پذیرفته و در آخر برای ارزیابی صحت 
لایه‌های ag‏ شده از دو روش دقت کل و ضریب LIS‏ استفاده شده است. 

2-۲-۲-۲ تصحیح هندسی و رادیومتریک 

در این مطالعه به‌منظور تصحیح هندسی» از نقشه‌های تویوگرافی با مقیاس 000( ag‏ شده از سازمان 
جغرافیایی ارتش استفاده گردید. در این مرحله با استفاده از نرم‌افزار 3 اظ اقدام به اعمال تصحیح‌های هندسی بر 
روی تصاویر شد و تصویر سنجنده OLI‏ سال ۱۳۹۸ با استفاده از روش تصویر- وکتور» زمین مرجع شد. برای این 
کار از ۳۹ نقطه کنترل زمینی با پراکنش مناسب و در تقاطع جاده‌ها؛ آپراهه‌ها و ... استفاده گردید. تا مدل ریاضی که 
برای پیدا کردن ضرایب مجهول در معادله مورداستفاده قرار می‌گیرد. حطای کمتری داشته باشد. برای تبدیل مختصات 
تصوير تصحیح شده به تصویر تصحیح نشده از تابع درجه اول استفاده گردید و برای نمونه‌گیری مجدد ارزش 
پیکسل‌های تصویر تصحیح نشده از روش نزدیک‌ترین همسابه استفاده و درنهایت سنجنده OLI‏ با حطای RMSE‏ 


1 Ghosh et al 


سال دوازدهم ارزیابی تأثیر تغییرات کاربری اراضی روی فرسایش و .... wa‏ 


به‌منظور تعیین تغییرات علاوه بر هم مختصات بودن تصاویر باید ابعاد پیکسل‌های آن‌ها نیز OLS‏ که در این تحقیق 
اندازه همه پیکسل‌ها ۲۸۵ متر است. 

در ضمن. تصحیح رادیومتریک زمانی انجام می‌گیرد که از تصاویر چند زمانه یعنی تصاویری که مربوط به فصول 
یا سال‌های مختلف و یا سنجنده‌های مختلف هستند. استفاده شود. دو نوع تصحیح رادیومتریک وجود دارد. تصحیح 
رادیومتریک مطلق و تصحیح رادیومتریک نسبی. روش تصحیح رادیومتریک مطلق نیازمند ورود داده‌های مربوط به 
خصوصیات اتمسفریک و کالیبراسیون سنجنده است. انجام این تصحیح در اغلب موارد و بهعصوص برای داده‌های 
قدیمی کار بسیار سختی است (دیو و همکاران؛ ۲۰۰۲). در مقابل تصحیحات رادیومتریک نسبی با هدف کاهش 
متغیرهای اتمسفریک و غیرمنتظره در ke‏ تصاویر چندزمانه صورت می گیرد. یکی از روش‌های تصحیح رادیومتریک 
نسبی» کاهش تیرگی پدیده" است. در حالت ایده‌آل پدیده‌های تیره‌رنگ دارای تابش صفر در همه طول‌موج‌ها هستند. 
در این روش فرض می‌شود که در هر باند از تصویر می‌توان پیکسل‌هایی CHL‏ که polis‏ آن‌ها صفر یا نزدیک به یک 
است (مثل (OT‏ به‌این‌ترتیب اثر اتمسفری تابش انحرافی به‌صورت یک مقدار ثابت به پیکسل‌ها در هر باند اضافه 
می‌گردد. به همین دلیل جهت حذف خطای رادیومتریک بایستی ارزش پیکسل‌های هر باند از حداقل by DN‏ به 
alg‏ کے رھ کا ری دیک ووی سات yeas Cul‏ کسه د رای از sip‏ په کار کر 
می‌شود (چاوز و مکنون" AME‏ در این مطالعه از این روش جهت تصحیح رادیومتریک تصاویر ذکر شده در بالا 


استفاده شده است. شکل ۲ مراحل مفهومی تحقیق را نمایش می‌دهد. 


1 Duet al 
2 Dark Object Subtraction 
3 Chavez and Mackinnon 


۳۸۰ جغرافیا و مخاطرات محبطی شمارة سوم 


استانداردسازی 


کاربری اراض ی /بوشش گیاهی 


و تصحیح bay‏ 
در ArcGIS‏ 


اجرای مدل EPM‏ برای دوره‌های مختلف 


نقشه فرسایش و رسوب برای دوره‌های مختلف 


شکل ۲- فلوچارت مراحل انجام کار 

۳-۲-۲- بررسی میزان فرسایش و رسوب با استفاده از مدل EPM‏ 
آگاهی از وضعیت فرسایش و تولید رسوب در حوضه‌های آبخیز نیاز به مطالعه و بررسی عوامل مؤثر در این 
dls‏ ده او کک Sits‏ ان کر ھا ر کی GS AO icles‏ نا گیری ,وس رب کر 
کشور چنین مطالعاتی با استفاده از روش‌ها و مدل‌هایی که در سایر کشورها توسعه یافتهاند. انجام می‌گیرد. در این 
تحقیق برای برآورد شدت فرسایش و تولید رسوب حوضه آبخیز موردمطالعه از روش EPM‏ و تکنیک‌های RS‏ و 
GIS‏ استفاده شده است. بعد از وارد کردن نقشه‌ها به محیط نرم‌افزار ARCGIS‏ و زمین مرجع کردن آن‌هاء لایه‌های 
مورد نیاز رقومی و پایگاه اطلاعاتی برای آن‌ها ایجاد شد.با استفاده از مدل EPM‏ عوامل موثر در فرسایش در 
زیرحوضه‌ها ارزیابی و امتیازبندی شد. در این روش که در حقیقت ضریب شدت فرسایش است. چهار عامل ضریب 
شدت فرسایش حوضه آبخیز CP)‏ ضریب کاربری اراضی (Ma)‏ ضریب حساسیت سنگ و خاک به فرسایش (Y)‏ و 
شیب متوسط حوضه آبخیز (1) در زیر حوضه‌های مختلف و درنهایت در کل حوضه آبخیز موردمطالعه قرار گرفت. 
در هر یک از این زیرحوضه‌ها بر اساس چهار عامل, مقدار 2 یا ضریب شدت فرسایش. متوسط سالانه رسوب ویژه 
(میزان رسوب تولید شده در واحد سطح) و مقدار رسوب تولید شده در طول سال در سطح زیرحوضه‌ها و بالاخره 


در کل حوضه محاسبه گردیده است. 
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به ۷ زیرحوضه تقسیم شد. سپس عوامل مؤثر در تولید رسوب در هر یک از زیرحوضه‌ها موردمطالعه قرار 
گرفت. در مرحله بعد میزان فرسایش و تولید رسوب با استفاده از فرمول‌ها و ضرایب مربوطه به دست آمد که 


۱-۳-۲-۲- تعیین ضریب استفاده از زمین یا ضریب کاربری زمین (XA)‏ 

بعد از انتخاب ترکیب باندهای مناسب برای هر دوره (11432, ۳717۷۲432 و 01543 برای تهیه نقشه کاربری 
اراضی نیاز به نمونه‌های تعلیمی از منطقه مطالعاتی است. نمونه‌های تعلیمی از طریق تفسیر بر روی تصاویر ماهواره‌ای 
و از روش ROI‏ در محیط نرم‌افزار ENVI‏ به دست می‌آید. نمونه‌های اخذ شده از تصاویر» Gyre‏ طبقات کاربری 
اراضی است که در فرایند پردازش تصاویر به‌عنوان پیکسل‌های آموزشی عمل می‌کنند. با توجه به اطلاعات 
جمع‌آوری شده و نمونه‌های آموزشی, نقشه‌های کاربری اراضی برای منطقه موردمطالعه به کمک طبقه‌بندی نظارت 
شده با روش بیشترین احتمال تهیه شده است. 

روش طبقه‌بندی حداکثر احتمال هنوز هم یکی از رایج‌ترین روش‌های طبقه‌بندی تصویر است (بولستد و لیسند " 
۱ در این روش طبقه‌بندی براساس واریانس و کوواریانس انجام می‌شود. در روش مذکور فرض بر این است که 
daa‏ مناطق آموزشی hl‏ پراکنش نرمال هستند. در حقیقت نمونه‌های کلاس‌های آموزشی باید معرف آن کلاس 
باشند بنابراین تا حد امکان باید از تعداد نمونه‌های بیشتری استفاده شود تا تغییرات بسیاری از ویژگی‌های طیفی در این 
گسترة پیوسته قرار گیرد. در طبقه‌بندی حداکثر احتمال» پیکسل موردنظر به کلاسی نسبت داده می‌شود که بیشترین 
احتمال تعلق پیکسل به آن کلاس وجود دارد (علوی پناه و احسانی» ۲۰۰۳). 

امتیاز ضریب استفاده از زمین بین ۰/۱ تا ۱ متغیر است. با استفاده از نقشه کاربری اراضی تولید شده برای سال‌های 
موردمطالعه و شرایط استفاده از زمین برای هر یک از کاربری‌هاء امتیازی اختصاص داده می‌شود و سپس با استفاده از 
مساحت اشغال شده توسط هر کاربری در هر زیرحوضه میزان متوسط وزنی ضریب استفاده از زمین در هر زیر 
حوضه محاسبه گردید. 

۲-۳-۲-۲- ضریب حساسیت سنگ و خاک نسبت به فرسایش CY)‏ 

امتیاز این عامل بین ۰/۲۵ تا ۲ متغیر است و بسته به مقاومت سنگ در مقابل فرسایش درجات مختلف به آن تعلق 
می گیرد. بدین ترتیب میزان امتیاز ضریب حساسیت خاک نسبت به فرسایش هر یک از سازندهای موجود در سطح 


حوضه از طریق جداول مربوطه تعیین گردید و سپس با توجه به وزن فرسایشی هر سازند و مساحت اشغال شده 


1 Region of Interest 
2 Bolstad and Lillesand 


YAY‏ جغرافیا و مخاطرات محبطی شمارة سوم 


توسط آن در هر زیرحوضه میزان متوسط وزنی ضریب حساسیت خاک نسبت به فرسایش در هر زیرحوضه محاسبه 
گردید. 

۳-۳-۲-۲- ضریب فرسایش حوضه CF)‏ 

امتیاز این عامل بین ۰/۱ تا ۱ متغیر است. این عامل در اصل وضعیت فعلی فرسایش را در سطح حوضه نشان 
می‌دهد. با استفاده از نقشه کاربری اراضی تولید شده برای سال‌های موردمطالعه و ویژگی‌های فرسایش هر کاربری» 
به هر یک از کاربری‌ها امتیازی اختصاص داده شده و سپس با استفاده از مساحت اشغال شده توسط هر کاربری در 
هر زیرحوضه. میزان متوسط وزنی ضریب فرسایش در هر زیرحوضه محاسبه شد. 

1-۳-۲-۲ تعیین شیب متوسط حوضه T‏ 

در این مطالعه جهت استخراج پارامتر شیب از نقشه مدل رقومی ارتفاع در محیط نرم‌افزار ۵0015 نقشه در 
صد متوسط شیب استخراج گردید. 

با توجه به چهار فاکتور محاسبه شله در مرحله قبل» مقدار 2 یا ضریب شدت فرسایش برای هر یک از 
زیرحوضه‌ها و درنهایت کل حوضه آبخیز با استفاده از daly‏ ۱ محاسبه گردید. 

Z=Xa.Y(¥+ 1%) \ رابطه‎ 

در روش EPM‏ برای تخمین متوسط سالانه فرسایش ویژه که در حقیقت میزان فرسایش را در واحد سطح 

(هکتار یا کیلومترمربع) حوضه را نشان می‌دهد. از رابطه Y‏ محاسبه می‌شود. 
رابطه ۲ Wsp = T.H.Z Pi‏ 

که در آن Wsp‏ مقدار فرسایش برحسب مترمکعب در سال در کیلومترمربع. H‏ متوسط بارندگی سالیانه حوضه 
برحسب میلی‌متر Z‏ شدت فرسایش. =Pi‏ عدد ثابت پی معادل ۲/۱۶ و T‏ ضریب درجه حرارت که از رابطه‌ی ۲ به 
دست می‌آید (احمدی ۲۰۰۷). 

T= (t/10 + 0.1( 5 Y رابطه‎ 

در این رابطه (t)‏ متوسط درجه حرارت سالیانه حوضه برحسب درجه سانتی گراد است. 

مقدار Wsp‏ محاسبه‌شده در رابطه ۲ بیان‌کننده مقدار SUE‏ است که از بستر خود جدا شده و انتقال یافته است» 
اما همه این مواد فرسایش Al‏ به حروجی حوضه نمی‌رسد و آن مقدار خاک فرسایش al‏ که به رسوب خروجی 
تبدیل خواهد شد. بستگی به متغیرهای مختلفی دارد که درواقع بر ضریب رسوب‌دهی حوضه تأثیر می‌گذارد (قنبرزاده 
و غلامرضایی» ۲۰۰۷). 

مقدار رسوب ویژه حوضه آبخیز از طریق رابطه 6 به دست می‌آید که در آن Gap‏ رسوب ویژه برحسب 


مترمکعب در سال در کیلومترمربع ۷ فرسایش ویژه برحسب مترمکعب در سال در کبلومترمربع Rus‏ ضریب 


سال دوازدهم ارزیابی تأثیر تغییرات کاربری اراضی روی فرسایش و .... YAY‏ 


که چه میزان مواد رسوبی از نقطه حروجی خارج می‌شود و مقدار آن از طریق رابطه ۵ حاصل می‌شود. در این رابطه 


می‌آید (رفیعی» (Yq‏ 


Gsp = We. Ru T 
Ra = 4P * (3 /L +10 

۵ رابطه‎ 
D=D,,-D 

رابطه \ 


در رابطه Dav A‏ ارتفاع متوسط حوضه آبخیز و Do‏ ارتفاع نقطه خروجی رودخانه است. مقدار رسوب کل 


حوضه از رابطه V‏ به دست می‌آید (احمدی. (Ye aV‏ 


= وت‎ . ۸ 
Gs sP ۷ رابطه‎ 


که در آن: Gs‏ رسوب کل برحسب مترمکعب در سال در کیلومترمربع؛ Gp‏ رسوب ویژه برحسب مترمکعب در 


سال در کیلومترمربع و ۸ مساحت حوضه آبخیز پرحسب کیلومترمربع است. 


۳- نتایج و بحث 
پس از agi‏ تصاویر ماهواره‌ای با استفاده از bla‏ کنترل زمینی» تصاویر با خحطایی حدود ۰/۳۲ زمین مرجع گردید 
تصاویر تطابق یافته. دقت بالای عمل تطابق هندسی مورد 
تائید قرار گرفت. همچنین برای تصحیح رادیومتریک از روش تصحیح رادیومتریک نسبی استفاده گردید. این روش 
برای کاهش اثرات پخش اتمسفریک بر روی تصویر است. در ضمن. تعداد پوشش/ کاربری اراضی با توجه به 
تصاویر و نقشه‌های کاربری موجود و شرایط منطقه موردمطالعه و بازدید میدانی برای تهیه نقشه پوشش زمین به 
چهار کلاس تقسیم شده است که عبارتند از: مرتع» کشاورزی؛ مسکونی و پهنه‌های آبی. 
پس از طبقه‌بندی کاربری‌های اراضی. به ارزیایی صحت پرداخته شد. برای ارزیابی صحت طبقه‌بندی, انتخاب 
یک‌سری پیکسل‌های نمونه معلوم و مقایسه کلاس آن‌ها با نتایج طبقه‌بندی لازم است. بدین‌جهت تصاویر طبقه‌بندی 
شده با نمونه‌های زمینی تهیه‌شده به نرم‌افزار ENVI‏ انتقال یافته و میزان دقت کاربری‌ها (دقت کل و ضریب کاپا) به 
دست آمد. دقت کل در نقشه‌های طبقه‌بندی شده سال‌های ۰۱۳۹۲ ۱۳۸۹ و ۱۳۹۸ به ترتیب ٩۵‏ ۹۸ و ۹۵ در صد 
بود. با توجه به نتایج حاصل از ارزیابی دقت در سه دوره موردمطالعه. میزان تطابق نمونه‌ها با واقعیت زمینی بسیار 
بالاست. 


TAF‏ جغرافیا و مخاطرات محیطی شمارة سوم 


۱-۳- انواع کاربری اراضی و میزان تغییرات آن‌ها 

در شکل‌های bY‏ ۵ نقشه‌های کاربری اراضی سال‌های ۱۳۹۹ ۱۳۸۹ و ۱۳۹۸ به ترتیب OLY‏ داده شده است. در 
ده ای celles‏ روط اا a A‏ هه اس هش hes‏ تیور 
درصد تغییرات طبقات پوشش زمین مرتع» کشاورزی» مسکونی و پهنه‌های آپی نشان داده شده است. چنانچه در 
شکل 1 مشاهده می‌شود. در تمام سال‌های موردمطالعه. بیشترین طبقه پوشش زمین مربوط به مرتع و کمترین آن 
مربوط به کاربری کشاورزی» مسکونی و پهنه‌های آبی است. همچنین در شکل 1 مشاهده می‌شود. در تمام دوره‌های 


موردمطالعه روند تغییرات مرتع Gy ped‏ کاهشی ولی روند تغیبرات کشاورزی به‌صورت افزایشی بوده است. 


O ke PAF. Pins papet, g °° YAT E ۲۳ 
| N Zo N < 
۹ L > BA g- 
. | S ۰ ۳ 
$ F 
۹ > g- 
I 7 کاربری اراضی‎ | eee 3 
| زیری اراضی | ج 2 رده‎ : 
55 ci: E کشاورزی‎ 
: M 91 ۳ 
SMe S> Bux > 
$- y 9%- g- 
6 3 0 6 6 3 0 6 
لس" "سس‎ i sls کار سح‎ 
PFA.“ Af. Y¥Pooes PE... PFA ۱۳۵2۰۰۰ Pe... ۱۳۶۰۰ 
۱۳۸۹ شکل ۳- نقشه کاربری اراضی در سال ۱۳۹۲ شکل - نقشه کاربری اراضی در سال‎ 
۲ PEA. AF... ۳۶-۰۰ ۳۶۶ 
3] ¥ 3 
100 لا‎ W a س‎ 
90 ۰ ۰ 
80 0 S : 
3 3 
70 > > 
60 ۹ = 
50 
40 : ۱ 1 
: - کاربری اراضی‎ : 
> > 
20 ۹ E کناررزی‎ pa 
۳ E £- 
۱ : : 
E è 
۹ ۹ 
بینه های آبی‎ 6 3 0 6 
| سس ۳۳ سس‎ | ie 
m1366 m1389 8 YFA Yaf ۳۶ Ee 
شکل ۵- نقشه کاربری اراضی در سال شکل *- نمودار تغییرات مساحت کاربری‎ 


۴۹۸ اراضی در بازه زمانی بین ۱۳۱ تا ۱۳۹۸ 


سال دوازدهم ارزیابی تأثیر تغییرات کاربری اراضی روی فرسایش و ... YAO‏ 
حدول ۱- مساحت کاربری‌های اراضی در سال‌های مختلف 
سال ۱۳۱۱ سال ۱۳۸۹ سال ۱۳۹۸ 
کاربری سال 
هکتار درصد هکتار درصد هکتار درصد 
کشاورزی ۱۳ 12 Y/Y‏ ۷/۳۱ ۱۳۵۹۹ ۳/۳۲ 
مرتع 101/0 0/4 0/۸1 VN\Y \Y¥ov/T\ E‏ 
مسکونی NA ۱۳۳/۳۳ ۰/۷۷ ۱۳۳/۳۰ ۰/۷۷ wy‏ 
پهنه آبی ۱/۸ LA ۱/۸ [A‏ 


مع‌الوصف بر اساس نتایج حاصل از این تحقیق» تغیبرات انواع کاربری اراضی طی بازه‌های زمانی ۱۳۲ و ۱۳۸۹ 
و ۱۳۹۸ به شرح زير است: 

۱-۱-۳- اراضی با کاربری کشاورزی 

سطح این اراضی در فاصله زمانی بین ۱۳۹۲ تا ۱۳۹۸ به میزان ۱7۷۸۳ درصد افزایش داشته است. بر اساس نتایج 
به‌دست‌آمده مشخص گردید سطح این اراضی ۷۱۷/۱۳ هکتار (۶/۱۹ درصد) در سال ۱۳۹۸ به ۱۲۳۲/۲۷ هکتار 
(۷/۲۱ درصد) در سال ۱۳۸۹ و ۳۵۹۳/۲۹ هکتار (۲۱/۰۲ درصد) در سال ۱۳۹۸ افزایش داشته است. از علل مهم 
افزایش وسعت این نوع کاربری می‌توان به پدیده مهاجرت و ازدیاد جمعیت در سال‌های اخیر و متعاقب آن تبدیل 


مراتع و اراضی BE‏ کاربری به باغات و اراضی کشاورزی اشاره نمود. 


2-۳-۳ اراضی با کاربری Er‏ 

سطح این اراضی در فاصله زمانی بین ۱۳۱۲ تا ۱۳۹۸ به میزان ۱۷۹۲ درصد کاهش داشته است. بر اساس نتایج 
به‌دست‌آمده مشخص گردید سطح این اراضی ۱۲۵۱/۳۵ هکتار (۹۵/۰۶ درصد) در سال ۱۳۹۲ به NOVY VAT‏ هکتار 
٩۱/۹۶(‏ درصد) در سال ۱۳۸۹ و ۱۳۳۵۷/۳۳۲ هکتار (۷۸/۱۲ درصد) در سال ۱۳۹۸ کاهش داشته است. دلیل این 
کاهش, تبدیل اراضی مرتعی به CEL‏ و اراضی کشاورزی در طی سالیان ۱۳۷۲ تا ۱۳۹۸ بوده است که خود از 


پیامدهای افزایش جمعیت در سال‌های اخیر است. 
۳-۱-۳- اراضی با کاربری مسکونی 
به‌دست‌آمده مشخص گردید سطح این اراضی ۱۳۱/۱۳ هکتار (۰/۷۷ درصد) در سال ۱۳۹۱ به ۱۳۲/۲۰ هکتار (۰/۷۷ 


YAF‏ جغرافیا و مخاطرات محیطی شمارة سوم 


1-۱-۳- اراضی با پهنه cal‏ 

سطح این اراضی در فاصله زمانی بین ۱۳۸۹ تا ۱۳۹۸ هیچ‌گونه تغییری نيافته است. نتیجه حاصل از بررسی 
نقشه‌های کاربری اراضی مربوط به بازه زمانی ۱۳۱۲ تا ۱۳۹۸ نشانگر تغییراتی در مساحت انواع کاربری‌ها است. 
ازجمله دلایل کاهش سطح مراتع می‌توان افزایش روزافزون جمعیت که موجب تغییر کاربری مرتعی و تبدیل آن به 
زمین‌های زراعتی» مناطق مسکونی و ... شده است. البته عوامل طبیعی دیگری نظیر عوامل اقلیمی مانند سرمای 
زودرس» خشک‌سالی و عوامل بیولوژیک مانند olab‏ آفات نیز در تخریب مراتع حوضه نقش بسزایی دارد؛ اما کاهش 
سطح مراتع حوضه را می‌توان به دلیل استفاده بی‌رویه دامداران جهت تعلیف احشام خود دانست؛ بنابراین با استفاده 
از تصاویر ماهواره‌ای می‌توان تا حدودی از وضعیت کاربری‌ها پیشین آگاه شد و این توانایی را به محقق می‌دهد که 
بتولند با تعیین کاربری‌ها و مقایسه با کاربری‌های فعلی تا لندازه‌ای از تغییرات مساحت‌ها و نقش این تغییرات در 
افزايش فرسایش و رسوب و حتی مشکلات عدیده‌ای نظیر سیل و... مطلع شد. 

۲-۳- مقادیر ضریب استفاده از زمین یا ضریب کاربری زمین (Xa)‏ 

مطالعه کاربری اراضی حوضه کال آجی OLE‏ می‌دهد. حوضه موردمطالعه از ٤‏ کاربری تشکیل‌شده است که در 
جدول ۲ انواع مختلف کاربری و در شکل‌های ۸۰۷ و ٩‏ ضرایب مربوط به هر نوع کاربری در سال‌های ۱۳۹۱ ۱۳۸۹ 
و ۱۳۹۸ و در جدول ۳ ضرایب مربوط به هر نوع کاربری در زیر حوضه‌های مختلف نشان داده شده است. کاربری 
اراضی به‌عنوان یکی از فاکتورهای اصلی فرسایش خاک محسوب می‌شود و با توجه به نقشه کاربری اراضی درصد 
خیلی زیادی از مساحت حوضه به کاربری مرتع و کشاورزی اختصاص داده شده است که امتباز حاصله در جدول 


ذیل نشان داده شده است. 


حدول ۲- انواع مختلف کاربری در حوضه و ضرایب استفاده از زمین 


شماره نوع کاربری ضریب استفاده از زمین (Ka)‏ 
۱ مسکونی ۰/۱ 
۲ کشاورزی 4/ 
مرتع ۵/ 
t‏ سطوح آبی ۰/۱ 


سال دوازدهم ارزیابی تأثیر تغییرات کاربری اراضی روی فرسایش و .... YAV‏ 
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ضریب استفاده از زمین (Xa)‏ 


" 0.9 
Low: 0.1 


۴۱۶۵۰۰۰ 
Fyd 


YFA ++ ۳۵۴ ۳۶۰ ۳۶۶۰ 


شکل -٩‏ نقشه ضریب استفاده از زمین در سال ۱۳۹۸ 


جدول ۳- مقادیر ضریب استفاده از زمین (Ka)‏ در هر یک از زیرحوضه‌های حوضه آبخیز کال آجی 


زیر حوضه‌های آبخیز ۱۳۵۰ ۳۵- \aaXa-‏ میزان تغییرات 
\ ۰/0۱ ۰/۲ ۰/۱ ۰/۱ 
y‏ ۰/۵۱ ۰/6۲ 00/* £" 
Y ۰۹ ۰/۹ ۳‏ 4/ 
3 ۰/0۱ ۰/0۱ ۰/5۹ ۵/ 
tA ۵‏ ۰:۹ \0/+ ۳/ 
st ۰/۲ tA ۸ 1‏ 
"0Y ۰0۰ ۰/1۹ ۷‏ £“ 


YAA‏ جغرافیا و مخاطرات محبطی شمارة سوم 


بر اساس جدول شماره ۳ می‌توان نتیجه گرفت مقدار ضریب استفاده از زمین در مجموع در همه زیرحوضه‌ها 
افزایشی بوده و در ضمن در زیرحوضه شماره ۱ مقدار آن حداقل (۰/۱) و در زیرحوضه شماره ۳ مقدار آن حداکثر 
)4/4( است. 

۳-۳- ضریب حساسیت سنگ و خاک به فرسایش CY)‏ 

امتیاز عامل زمین‌شناسی از منابع مختلف (تنگستانی» VEN‏ یمانی و اسکندرنژاد. ۲۰۱۰؛ خاکسار و همکاران؛ » 
۰۱ اخذ و طی آن از نظر مقاومت برای سنگ‌شناسی دامنه بین ۰/۱ تا ۲ BES‏ شده است. به‌طوری‌که سنگ‌های 
سست و بسیار حساس به فرسایش بیشترین امتیاز (کد ۲) و سنگ‌های مقاوم و سخت» کم‌ترین امتیاز AS)‏ ۰/۱) را به 
خود اختصاص داده‌اند. با توجه به گزارش زمین‌شناسیء منطقه موردمطالعه از لحاظ لیتولوژی و حساسیت به فرسایش 
به ۷ کلاس حساسیت تفکیک گردیده که در شکل ۱۰ مقادیر ضریب حساسیت سنگ و خاک به فرسایش در حوضه 
آبخیز IIS‏ آورده شده است. 

۳-- محاسبه شیب متوسط D‏ 


شیب حوضه رابطه مستقیم و نسبتاً پیچیده‌ای با نفوذ آب و جریان سطحی و نیز رطوبت خاک دارد و در میزان 
فرسایش و تولید رسوب حوضه موثر است؛ بنابراین وجود این عامل به‌صورت بحرانی نقش تأثیرگذار در روند 
فرسایش و تولید رسوب می‌تواند داشته باشد bs)‏ و لی ۲۰۰۲). دامنه کلاس‌های شیب به‌کاررفته در این مطالعه 
برحسب درصد عبارتند از: ۵-۰ ۱۲-۵ ۲۵-۱۲ و ۰-۲۵. در شکل ۱۱ نقشه شیب منطقه و همچنین در جدول ٤‏ 
شیب حوضه آبخیز کال‌اجی نشان داده شده است. همان‌طور که از جدول ۶ پیداست بخش عمده حوضه موردمطالعه 
دارای شیب بیش از ۱۲ درصد VO)‏ درصد مساحت حوضه) داشته‌اند که نشان از اهموار بودن حوضه و تأثیر ol‏ 


1 Dai and Lee 


YAN 


سال دوازدهم اززیابی تأثیر تخییرات کاربری اراضی روی فرسایش weg‏ 
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شکل ۱۱- نقشه درصد شیب حوضه 
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شکل ۱۰- نقشه ضریب حساسیت سنگ و 
خاک به فرسایش (Y)‏ 


جدول .٤‏ شیب حوضه 


طبقات شیب (درصد) مساحت (هکتار) درصد مساحت 
۵-۰ 2۳۵/۹ ۱/۹۶ 
WAY/V ۱۲-۵‏ ۳۹/۳۲ 
۲۵-۲ 0441/14 ۳۵/۰۳ 
£+-Yo‏ 0/۷1 ۰/۹ 


poles -۵-۳‏ ضریب فرسایش منطقه )¥( 


پس از انجام بازدیدهای متعدد از حوضه موردمطالعه» نقشه اشکال مختلف فرسایش gba wes Fangs‏ رکلی در 


حوضه کال‌آجی ۵ تیپ فرسایش مشاهده گردید که در جدول ۵ انواع مختلف فرسایش و ضرایب مربوطه آمده 


است. همچنین در شکل 11 به‌صورت نقشه ضریب فرسایش حوضه نشان داده شده است. 


جدول ۵-انواع مختلف فرسایش در حوضه و ضرایب فرسایش 


۲ شماره نوع فرسایش ضریب فرسايش‎ 
۷ Dulin ۱ 
vas River Bank Erosion ۲ 


۳۹۰ جغرافیا و مخاطرات محیطی شمارة سوم 


Ve Gully ۳ 
۱/۷ River Erosion + Sheet Erosion + Micro terrace £ 
1 sheet + rill + river + splash erosion ° 


24 N : 
> < 
+ |y E 
۰ 5 ۰ 
3] ۳2 
= Z 
سو س‎ 
> 
سو‎ 3- 

(psi) ضریب فرسایش حوضه‎ 
. High: 1 2 
i 6 E 
w Low : 0.3 ` 
g- سل‎ 

5 2.5 0 

کیلومتر 
Y¥A- ++ YO¥- ++ Ee ۳۶۶۰۰۰‏ 


شکل ۱۲- نقشه ضریب فرسایش حوضه 


uy pe"‏ شلات فوسایش 
پس از تهیه عوامل چهارگانه مدل و جای‌گذاری آن‌ها ضریب شدت فرسایش در اشکال AY‏ ۱۶ و ۱۵ در 
سال‌های ۰۱۳۱۱ ۱۳۸۹ و ۱۳۹۸ نشان osla‏ شده است. همچنین در جدول 7 شدت فرسایش در زیرحوضه‌های 
مختلف را می‌توان مشاهده کرد. Cle‏ شدید بودن فرسایش در این زیرحوضه‌ها» شیب زیاد و حساسیت سنگ‌ها 


نسبت به فرسایش و خصوصاً تغیبرات کاربری اراضی است. 


Aeee ۱۳۰۰ Fates EE... YFAcee Afe VF owes VFS ees 
1 N l mt: if RN 
2 کک‎ 
+ Iw. E + * | E 
. 5 Jf, S 
$ $ $] 
۹ g سو‎ 
: i 
9 g > 
8 ضریب شدت فرسایش(2)‎ 
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۲ وی‎ E High : 1.94068 
۰ High : 2.02222 ۰ E 
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+ £ + 
5 25 0 5 25 5 
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۱۳۸۹ سال‎ ۱۳۹ 
. YFA- YaF- ۳۶۰۰ ۳۶۶۰۰ 
a, N : 
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+ ly E 3- 
i 5 ۲ 
4 | 4 
A = 
و‎ + 
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4- g- 
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کیلومتر‎ 
۱۳۴۸۰۰ ۳۵۰۰۰ ا‎ EE... 


شکل ۱۵- نقشه ضریب شدت فرسایش (Z)‏ در سال ۱۳۹۸ 


۴۶۵۰۰۰ 


Yay‏ جغرافیا و مخاطرات محیطی شمارة سوم 


جدول *- ضریب شدت فرسایش زیر حوضه‌های آبخیز کال آجی 


زیر حوضه‌های آبخیز 7-۲ 7-۹ 7-۸ شدت فرسایش 
“AY ۰/۸۰ ۱‏ ۰۹۵ شدید 
AY AY ۲‏ ۰/۸۸ شدید 
“VA ۳‏ ۸۱ ۰/۹۲ شدید 
AY ۸۱ t‏ ۰/۸۹ شدید 
NY ۵‏ ۱/۷۳ ۳ شدید 
IN +/VO Né 4‏ شدید 
IN IN /۸ ۷‏ شدید 
میانگین “VA‏ ۸ ۰/۸۷ 


بر اساس جدول شماره ٦‏ می‌توان نتیجه گرفت مقدار ضریب شدت فرسایش حوضه در فاصله زمانی ۱۳۹۲ تا 
۸ در مجموع در همه زیرحوضه‌ها افزایشی بود و در ضمن در زیرحوضه شماره ۳ مقدار آن حداقل (۰/۱) و در 
زیرحوضه شماره ٤‏ مقدار آن حداکثر (۰/۸) است. 

۷-۳- برآورد فرسایش ویژه در هر یک از زیر حوضه‌ها 

در ادامه نتایج محاسبات برآورد مقادیر فرسایش ویژه هر یک از زیر حوضه‌ها در اشکال AT‏ ۱۷ و ۱۸ همچنین 
در جدول ۷ و ۸ و شکل ۱٩‏ درصد فرسایش ویژه هر یک از زیر حوضه‌های حوضه آبخیز کال‌آجی ارائه گردیده 
است. با توجه به جداول ۸ و ٩‏ فرسایش ویژه برحسب مترمکعب در سال در کیلومترمربع در همه زیرحوضه‌ها در 
فاصله زمانی ۱۳۱۲ تا ۱۳۹۸ کماکان افزایشی بوده منتها در زیرحوضه شماره ۵ مقدار آن حداقل AVON)‏ و در 
زیرحوضه شماره ۳ (۳۲۸/۹۸) مقدار آن حداکثر است. همچنین» نتایج حاصل از محاسبات نشان می‌دهد مقدار 
فرسایش ویژه هر یک از زیرحوضه‌ها در دوره‌های زمانی مختلف متفاوت است. به‌طوری‌که زیر حوضه شماره ۵ و ۶ 
در هر سه دوره نسبت به ple‏ حوضه‌های دیگر مقدار فرسایش بیشتری را به خود اختصاص داده است و زیرحوضه 


وا چ که شماره ۱ gs‏ از وهای بیج ۱۳/۵ E‏ ۱۳۹۸ کین مقار فرسایقن را کرو اتضاصن استم 


سال دوازدهم ارزیابی تأثیر تغییرات کاربری اراضی روی فرسایش و vay ma‏ 


° Maece ۳۰۰ r eta ۳۴۳۰۰۰ FA PAS.“ aE 
= | 5I > ۱ f L 3 
z| À ا‎ | š 
£ i £ $ ۹ 
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z) } 3 
$ + 
W E 
: 5 : 
متوسط سالانه فرسایش ویژه > و‎ E 
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شکل ۱۸- نقشه فرسایش ویژه هر یک از زیرحوضه‌ها در سال ۱۳۹۸ 


حدول -V‏ درصد فرسایش ویژه هر یک از زیر حوضه‌ها در بین بازه زمانی ۲ تا ۱۳۹۸ 


Wsp (m/ha/yr) T : 
زیر حوضه‌های ابخیز‎ 
۱۳۹۸ ۱۳۸۹ ۱۳۹ 
0۸ ۱۳/۳ ۱۳/۳ ۱ 
۱۳/۳۱ ۱۳/۹۶ ۱۳/۹۶ ۲ 
۱۳/۳۱ ۱/۳۶ ۱۱/۳۹۶ ۳ 


۳۹۴ جغرافیا و مخاطرات محیطی شمارة سوم 


Wsp (m*/ha/yr) 5 ۱‏ 
زیر حوضه‌های ابخیز 

۱۳۹4۸ ۱۳۹ ۱۳۹۹ 

۱۳/۵۰ ۱۳/۳۰ ۱۳/۳۰ 1 

۱/۸۶ ۱/۸ ۱۰/2۸ ۵ 

۱/0۸ ۱۱/۳۹ ۱/۳۹ 1 

۱۳/۸۲ ۱۳/۹ ۱۳/۹۶ ۷ 

میانگین ۱۳/۳ ۱۳/۶۲ ۱۳/۹۹ 


جدول ۸- فرسایش ویژه هر یک از زیر حوضه‌های حوضه آبخیز کال آجی در بین بازه زمانی ۱۳۹۳ تا ۱۳۹۸ 


۱۷۷۲5۵۰۱۳۹۸ (M?/Km7/Year) ۱۱۷۷۲5۵-۱۳۸۹ (MP/Km7/Year) ۱۷۷5۵۰۱۳۰۱۰ (M3/Km7/Y ear) زیر حوضه‌های آبخیز‎ 

V/A: \YVA/YV 10/۹۰ \ 

\EEV VE ۱۳۳۰/۵ ۱۳/۳ ۲ 
\EATIEY 1/۰97 ۱۱۱۶ ۳ 
00/۸1 TEVE ۱۳۳/۳۵ 3 
WUT VAYA ۱۰/۹ ۵ 
AKATAA 107/4 A4۲ 1 
\ore/yy \YAY/\V (2 ۷ 
۱۰۰۹/۵ "1+۵ ۸/4 جمع‎ 
۱۶۳/۱ ۱۳9۹/۸ ۱۳۳۹/۰۹ میانگین‎ 


Wsp (m3/ha/r) 


٩ ۲ ۳ م‎ ò و‎ v 


m ۱۳۶۶ E ۱۳۸۹ E ۸ 


شکل -۱٩‏ نمودار درصد فرسایش ویژه هر یک از زیرحوضه‌های کال آجی در بین بازه زمانی ۱۳۹۲ تا ۱۳۹۸ 


سال دوازدهم ارزیابی تأثیر تغییرات کاربری اراضی روی فرسایش و ... 1۹۵ 


AT‏ برآورد رسوب ویژه در هر یک از زیرحوضه‌ها 

با توجه به مقادیر متوسط ضریب شدت فرسایش (Z)‏ در هر یک از واحدهای کاری حوضه آبخیز کالآجی» 
رسوب ویژه مربوط به هر زیرحوضه با استفاده از روش میانگین وزنی محاسبه گردید. در ادامه نتایج محاسبات برآورد 
polis‏ رسوب هر یک از زیر حوضه‌ها در اشکال ۰۲۰ ۲۱ و ۲۲ همچنین در جدول ٩‏ و ۱۰ و شکل ۲۳ درصد 
رسوب ویژه هر یک از زیرحوضه‌های حوضه آبخیز کال‌آجی ارائه گردیده است. با توجه به جدول ۱۰ رسوب ویژه 
برحسب مترمکعب در سال در کیلومترمربع در همه زیرحوضه‌ها به‌جز زیر حوضه شماره ۷ در فاصله زمانی ۱۳۹۸ تا 
۸ کماکان افزایشی بوده است. همچنین» geld‏ حاصل از محاسبات نشان می‌دهد مقدار رسوب هر یک از 
زیرحوضه‌ها در دوره‌های زمانی مختلف متفاوت است. به‌طوری‌که زیرحوضه شماره ۷ و ۱ در هر سه دوره نسبت به 
سایر زیرحوضه‌های So‏ مقدار رسوب بیشتری را تولید کرده است و زیرحوضه ۵ و زیرحوضه شماره 1 در بازه 


زمانی بین ۱۳۱۱ تا ۱۳۹۸ کمترین مقدار رسوب را تولید کرده است. 


"ekki mens ll, We ۱۳۳۹۰ var ۳۶۳ ۱ 
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۰ ۱ 6.84 ne ۱01 ۳ 
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£ 2: و لا‎ LEJ 1 
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Yese ۱۳۹۹۰۶ Fee es EYe. Yee PAS... EF... 
رسوب ویژه هر یک از شکا ۱- نقشه رسوب ویژه هر یک از‎ 23 ٣ شکا‎ 


زیرحوضه‌ها در سال ۱۳۹۲ زیرحوضه‌ها در سال ۱۳۸۹ 


var 


FlAceee 


PAYE. ° 


PEN. 


متوسط سالانه رسوب ویژه 
(مترمکعب درهکتار درسال) 


2.24 ۱ 
7.45 
[yy 7-62‏ 
7.84 ۲ 
تن ا 
6.69 ۲ 
7s‏ ا 


۳۵۶ 


ae 


FF Aces ۴۱/۴۰۰۰ ۴۱۸۰ 


EY... 


۳۴۳۳۰۰۰ 


۳۶۳ 


شکل ۲۲- نقشه رسوب ویژه هر یک از زیرحوضه‌های حوضه آبخیز کال آجی در سال ۱۳۹۸ 


جدول -A‏ رسوب ویژه هر یک از زیرحوضه‌های حوضه آبخیز کال آجی در بین بازه زمانی ۱۳۹۲ تا ۱۳۹۸ 


Gsp-\v4, (M/Kn7/Year) 6280-۱۳۸۵ (M?/Km?/Y ear) Gsp- iri (M?/Km7/Year) زیر حوضه‌های آبخیز‎ 
AYE/\E VOVAN Wo/+4 \ 
۱/۹۲ T/A: Wiser ¥ 
VIYA Wy V/A ۳ 
۱۷/۸/۹ ۳/۹ 0۹/71 t 
OW to OOTAV 2۱۲ ۵ 
Ywa. 040/۲ 0/44 1 
VAVA 71/۲0 ۸۸ ۷ 
0777/4۷ C00۰۰ YoY Y جمع‎ 


حدول 7 درصد رسوب ویژه هر یک از زیرحوضه‌ها در بین بازه زمانی ۲ تا ۱۳۹۸ 


Gsp (m*/ha/yr) : 
۱۳۹۸ ۱۳۸۹ ۱۳2 زیر حوضه‌های ابخیز‎ 
NYE voy VYO \ 
V/£0 VAL A ۲ 
V/V VYY VW ۳ 
۷/۷۳۶ V۹ 0/۹ t 
0/۹۸ 0/0۷ ZANA ۵ 
14 0/40 5/۱۳۰ 1 
۷/۸۱ VAL Na ۷ 


سال دوازدهم اززیابی تأثیر تخییرات کاربری اراضی رو فرسایش weg‏ ۳۹۷ 


Gsp (m3/ha/yr) 


\ ۴ hg ¥ ۵ F ۷ 


Eres Myrna E \yrqa 


شکل ۲۳- نمودار درصد رسوب ویژه هر یک از زیرحوضه‌های حوضه آبخیز کال آجی در بین بازه زمانی ۱۳۱۲ 
تا ۱۳۹۸ 


پر اساس نتایج به‌دست‌آمده مجموع مقادیر رسوب ویژه مربوط به بازه زمانی ۱۳۹۲ (مترمکعب 4۳۵۳/۱۲) بوده 
که این میزان در سال ۱۳۹۸ به ۵۱72/۰۰ مترمکعب افزایش داشته است که با تغییرات افزایشی یا کاهشی 
تغییرات کاربری اراضی رخ داده در منطقه همخوانی دارد. چنانچه مشاهده شد که سطح اراضی کشاورزی در 
فاصله زمانی بین ۱۳۹۶ تا ۱۳۹۸ به میزان ۱0/۸۳ درصد افزایش داشته است. 

با توجه به جدول ۱۰ و با در نظر گرفتن وزن مخصوص رسوب معادل ۱/۳ تن در مترمکعب در کل حوضه به 
ترتیب برای سال‌های ۱۳۹۲ ۱۳۸۹ و ۱۳۹۸ آورده شده است که مقدار رسوب برای بازه زمانی ۱۳۹۲ معادل 
۲ تن در سال» برای بازه زمانی ۱۳۸۹ ales‏ ۵۸۸۷/۷۲ تن در سال و برای بازه زمانی ۱۳۹۸ 
معادل 1۷۱۳/۲۰ تن در سال است. که با تقسیم بر مساحت حوضه مقداری برابر ۰/۳۳ تن در هکتار در سال 
۲ ۶ تن در هکتار برای سال ۱۳۸۹ و ۰/۳۹ تن در هکتار برای سال ۱۳۹۸ برآورد گردید. 


t‏ نتیجه گیری 

رسوب صورت گرفته است. نتایج طبقه‌بندی تصاویر در این بررسی So Whe‏ تحقیقات صورت گرفته 
SL‏ قابلیت تصاویر ماهواره‌ای در تهیه نقشه کاربری اراضی و پوشش زمین است. استفاده از تصاویر چند 
زمانه با بازه‌های زمانی معین» می‌تواند به مدیران در پایش منابع طبیعی و تصمیم‌گیری درباره آینده کمک LS‏ 


نقشه‌های کاربری اراضی تهیه شده از حورضه موردبررسی Olas‏ می‌دهد که درصد مساحت کشاورزی» ‘ey‏ 


YAA‏ جغرافیا و مخاطرات محبطی شمارة سوم 


مسکونی و پهنه‌های آبی در سال ATVI‏ به ترتیب EAA‏ ۵/۰۶ ۰/۷۷ و " بوده که در سال ۱۳۸۹ به ترتیب به 
۱ ۸۱/۹۶ ۰/۷۷ و ۱/۰۸ رسیده و در سال ۱۳۹۸ به ترتیب ۵۲۱/۰۲ ۸۷۸/۱۲ ۰/۷۸ و ۰/۰۸ رسیده است؛ 
به‌عبارت‌دیگر در ۳۲ سال گذشته (۱۳۹۸-۱۳۲۲) حدود ۲۸۷۲ هکتار از سطح مراتع کاسته شده که تقریباً تمام 
آن به کشاورزی bas‏ شده است. درصد کاهش مساحت مرتع در دوره زمانی ۱۳۸۹-۱۳۹ و ۱۳۹۸-۱۳۸۹ 
به ترتیب ۳/۱۰ و ۱۳/۸۲ است. سرعت رشد مسکونی در حوضه موردمطالعه بسیار کم است. بررسی تغییرات 
کاربری اراضی و پوشش زمین سال‌های مختلف در حوضه نشان داد که تغییرات صورت گرفته عموماً به‌صورت 
کاهش سطح مرتع و افزایش مساحت اراضی کشاورزی است. بر این اساس می‌توان نتیجه‌گیری کرد که عامل 
اصلی تغییرات کاربری اراضی در حوضه موردمطالعه. دخالت‌های انسانی است. می‌توان گفت نتایج این مطالعه 
با نتایج مطالعه افرادی چون اسفندیاری و همکاران (۲۰۱۶) مرادی و همکاران TOV)‏ فیضی‌زاده (۲۰۱۷ 
اصغری‌سراسکانرود و همکاران (۲۰۱۷). مارتینز و همکاران (۲۰۱۱) و داسیلوا و همکاران (YONA)‏ تطابق دارد. 
بدین‌صورت که آن‌ها نیز معتقدند کاربری کشاورزی از بیشترین پتانسیل و مراتع از کمترین پتانسیل برای ایجاد 
فرسایش برخوردارند. نتایج نشان داد که ضریب شدت فرسایش در سال‌های AYTI‏ ۱۳۸۹ و ۱۳۸۲ به ترتیب 
۸ ۰/۸۰ و ۰/۸۷ است. مشاهده می‌شود که روند تغییرات ضریب شدت فرسایش افزایشی است. در ضمن» 
میزان فرسایش ویژه حوضه در این سال‌ها به ترتیب ۱۲۲۵/۰۹ ۱۲۵۹/۰۸ و ۱۶۳۵/۸۱ مترمکعب در کیلومترمربع 
در سال است؛ یعنی سال ۱۳۹۸ متناظر با ضریب شدت فرسایش» فرسایش ویژه بیشتری داشته است: کشت کار 
روی سازندهای سیلتی-رسی حساس به فرسایش آبی» احداث تاکستان‌های دیم روی دامنه‌های شیب‌دار و 
تخریب مراتع منجر به فرسایش تشدیدی در خاک‌های سطح حوضه شده است. این تغییرات نقش مهم مراتع 
و پوشش گیاهی حوضه آبخیز کال آجی را نشان می‌دهد که موجب افزایش میزان رسوب خروجی زیر حوضه‌ها 
گردیده است. 

به نظر می‌رسد لازم است در جهت حفظ عرصه‌های طبیعی, تثبیت و قانونی کردن کاربری اراضی و کنترل 
فرسایش. اقداماتی برای حفاظت از خاک و آب در مناطق با پتانسیل فرسایش بالا در چهارچوب طرح‌های 
حفاظتی در دستور کار متخصصان و مدیران اراضی قرار گیرد. از سوی دیگر» هرگونه تغییر در کاربری اراضی 
حوضه موردمطالعه. بایستی بر اساس آگاهی و دانش کافی از محدودیت‌های زیست‌محیطی منطقه و خطرات 
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